UNCLASSIFIED 
AD  NUMBER 


AD525863 

CLASSIFICATION  CHANGES 

TO: 

UNCLASSIFIED 

FROM: 

SECRET 

LIMITATION  CHANGES 

TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors ; 

Administrative/Operational  Use;  NOV  1972.  Other 
requests  shall  be  referred  to  Rome  Air 
Development  Center,  Griffiss  AFB,  NY  13441- 
5700. 


_ AUTHORITY _ 

RADC  ltr  dtd  16  Jul  1974;  RADC  ltr  dtd  16  Jul 
1974 


THIS  PAGE  IS  UNCLASSIFIED 


SECURITY 

MARKING 


The  classified  or  limited  status  of  this  report  applies 
to  each  page,  unless  otherwise  marked. 

Separate  page  printouts  MUST  be  marked  accordingly. 


THIS  DOCUMENT  CONTAINS  INFORMATION  AFFECTING  THE  NATIONAL  DEFENSE  OF 
THE  UNITED  STATES  WITHIN  THE  MEANING  OF  THE  ESPIONAGE  LAWS,  TITLE  18, 
U.S.C.,  SECTIONS  793  AND  794.  THE  fRANSMISSION  OR  THE  REVELATION  OF 
ITS  CONTENTS  IN  ANY  MANNER  TO  AN  UNAUTHORIZED  PERSON  IS  PROHIBITED  BY 
LAW. 


NOTICE:  When  government  or  other  drawings,  specifications  or  other 
data  are  used  for  any  purpose  other  than  in  connection  with  a  defi¬ 
nitely  related  government  procurement  operation,  the  U.  S.  Government 
thereby  incurs  no  responsibility,  nor  any  obligation  whatsoever;  and 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any 
way  supplied  the  said  drawings,  specifications,  or  other  data  is  not 
to  be  regarded  by  implication  or  otherwise  as  in  any  manner  licensing 
the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use. or  sell  any  patented  invention  that 
may  ip  any  way  be  related  thereto. 


BEST 

AVAILABLE  COPY 


SECRET 


RADC-TR-73-94 
Technical  Report 
November  1972 


EXPERIMENTS  AND  MODELS  FOR  PRAIRIE  SMOKE  (U) 
Aeronomy  Corporation 


lO 


Sponsored  by 

Defense  Advanced  Research  Projects  Agency 
ARPA  Order  No.  1423 


hibited  by  law. 


a  views  and  conclusions  contained  in  this  docuaent  are  those  of  the 
thors  and  should  not  be  interpreted  as  necessarily  representing  the 
official  policies,  either  expressed  or  implied,  of  the  Defense  Ad¬ 
vanced  Research  Projects  Agency  or  the  U.  S.  Government. 


1 


I 


V 


e  ct  L v*'1*3 


Rome  Air  Development  Center 
Air  Force  Systems  Command 
Griffiss  Air  Force  Base,  New  York 


.t  'to 


SECRET 


OOC  CONTROf 
Z12G2 


Do  not  return  this  copy, 
security  regulations. 


When  not  needed,  destroy  in  accordwce  with  pertinent 


SECRET 


EXPERIMENTS  AND  MODELS  FOR  PRAIRIE  SMOKE  (U) 

S.  A.  Bowhlll 
E.  E.  Mendenhall 


Contractor:  Aeronomy  Corporation 
Contract  Number:  P30602-72-C-0214 
Effective  Date  of  Contract:  1  September  1971 
Contract  Explratlcn  Date:  1  October  1973 

.'mount  of  Contract:  $216,000.00 
Program  Code  Number:  1E20 


Principal  Investigator:  Dr.  S.  A.  Bowhlll 


Phone: 

217  359-8007 

Project  Engineer: 

Vincent  J.  Coyne 

Phone: 

315  330-3141 

Contract  Engineer: 

Richard  U.  Carman 

Phone: 

315  330-3144 

This  research  was  supported  by  the 
Defense  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense 
and  was  monitored  by  Richard  W.  Carman 
RADC  (OCSE) ,  GAFB,  NY  13441  under  Con¬ 
tract  F30602- 7 2-C-O  2 14 . 


Classified  by  Ivory  Coral  Security  Guide 
EXEMPT  FROM  GENERAL  DECLASSIFICATION 
SCHEDULE  OF  EXECUTIVE  ORDER  11652 
EXEMPTION  CATEGORY:  3 
DECLASSIFY  OR  31  December  1983. 

Copy  /*/-  of  14  copies 
#231 


SECRET 


DDC  CONTROL 
NO.  31262 


yam 


PUBLICATION  REVIEW 


This  technical  report  has  been  reviewed  and  is  approved. 


/ 


< 


PADC  Project  Engineer 


■'ropjr 


9 


_ .  [ 


▲ 


SECRET 


ABSTRACT  (U) 

This  report  describes  satellite  transmission  experiments  carried  out  in 
Prairie  Smoke  lb  and  II  to  investigate  the  horizontal  extent,  altitude  dis¬ 
tribution,  intensity  and  frequency  dependence  of  field-aligned  artificial 
spread  F  (ASF)  irregularities.  Using  these  data,  a  quantitative  morphological 
model  is  described  in  which  the  irregularities  are  contained  within  a  tilted 
cylinder,  with  axis  parallel  to  the  magnetic  field  direction,  and  a  radius  of 
50  km,  centered  30  km  north  of  the  heater  at  an  altitude  of  300  km.  The 
fractional  fluctuation  of  electron  density  was  found  to  be  about  0.7  percent 
at  full  power,  and  the  correlation  distance  normal  to  the  field  was  about 
100  m.  The  irregularities  were  found  to  extend  through  the  entire  F  region 
between  200  and  400  km  altitude.  The  depth  was  found  to  vary  inversely 
with  frequency,  in  accordance  with  theory.  Alto  presented  in  this  report 
is  a  model  for  the  small-scale  structure,  including  an  analytic  technique  to 
determine  the  thickness  of  the  disturbed  region  from  the  change  in  echo 
return  range  with  elevation  of  a  monostatic  radar;  semi-thicknesscs  of  about 
25  km  are  found. 
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1.  INTRODUCTION  (U) 

(S)  Initial  experiments  on  satellite  signals  transmitted  through  a  dis¬ 
turbed  ionospheric  region  over  the  heater  at  Platteville,  Colorado,  were 
described  in  a  previous  report  (Bowh  '  and  Mendenhall,  1972).  The  present 
report  describes  the  continued  program  of  measurements  and  interpretation, 
primarily  in  Prairie  Smoke  lb  (February  1972)  and  II  (May  1972).  These 
results  are  presented  and  discussed  in  Sections  2  through  6. 

(S)  During  the  summer  of  1972,  a  scattering  model  for  the  artificia  -P 
F  (ASF)  irregularities  was  developed,  and  is  presented  in  Sections  7  and  1 1 . 

In  addition,  a  simple  model  for  on-frequency  scatter  (OFS)  has  been  developed 
and  is  presented  and  evaluated  in  Sections  8  through  10. 

(S)  In  preparation  for  possible  applications  studies,  some  considerations  for 
bistatic  OFS  modeling  are  given  in  Section  12. 
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2.  OBJECTIVES  OF  PRAIRIE  SMOKE  lb  AND  II  (U) 

(S)  Two  type?  of  satellite  transmission  experiments  have  been  carried  out 
by  Aeronomy  Corporation  in  connection  with  the  Prairie  Smoke  activities;  a 
geostationary  satellite  and  an  orbital  satellite  experiment.  The  ^wo  experi¬ 
ments  provide  a  means  for  determining  a  morphological  model  (and  hence,  a 
radar  scattering  model)  and  for  ultimately  determining  a  yield  model  of  arti¬ 
ficial  spread-F  (ASF)  irregularities  and  on-frequency  scatter  (OFS)  produced 
by  the  ionospheric  heating. 

(S)  A  detailed  list  of  objectives  for  the  satellite  experiments  was  pre¬ 
sented  by  Bowhill  and  Mendenhall  (1972).  Table  2.1  illustrates  the  appli¬ 
cability  of  the  results  of  experiments  conducted  during  Prairie  Smoke  lb  and 
II  for  satisfying  these  objectives.  Sections  2.1  and  2.2  describe  the  specific 
objectives  of  Prairie  Smoke  lb  and  II  experiments. 

2.1  Objectives  of  Prairie  Smoke  lb  (U) 

(S)  Due  to  the  relatively  short  time  available  to  mount  a  field  experiment 
for  Prairie  Smoke  I,  attention  was  concentrated  on  implementing  a  mobile  geo¬ 
stationary  satellite  receiving  station.  An  orbital  experiment  was  conducted 
on  December  6,  1971  when  a  single  pass  was  recorded  from  a  fixed  site  in  Pine 
Bluffs,  Wyoming.  Prairie  Smoke  lb,  which  was  scheduled  for  February  14-17, 
1972,  however,  offered  the  first  opportunity  to  deploy  a  mobile  facility  for 
recording  appropriate  orbital  passes.  Because  of  the  remoteness  of  the  site 
location  for  the  geostationary  experiment  and  the  potentially  inclement 
weather,  attention  was  concentrated  on  the  orbital  experiment. 

(S)  The  purposes  for  going  into  the  field  at  this  time  were  to  test  hard- 
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Table  2.1 


Potential  (P)  or  demonstrated  (D)  applicability 
tcllite  transmission  experiment  to  the  various 
using  amplitude  measurements  from  1-3  antennas, 
teroferometric  pair  (I). 


of  the  sa- 
objectives 
or  an  in- 


Ob  -j  ective 


Geostationary 


Orbital 


1  pi  eld  alignment 
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3  I 
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?  Size  and  orientation 
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D 
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3  Altitude  distribution 

D 

D 

4  Horizontal  extent 
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5  Onset  and  decay 

D 

D 

D 

A  Drift  velocity 
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7  Padinp  with  frequency 
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«  Ana l p,  flutter  with  frequency 

P 

9  Multipath  effects 

P 
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in  -pT*nni  mnCV  scatter 
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11  Modulation  degradation 

P 

P  _ 

n  Upflfpr  n^lrameters 
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i x  of f -vertical  heating 
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14.  New  antenna  concepts 
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ware  planned  for  future  orbital  experiments,  to  demonstrate  the  usefulness 
of  orbital  measurements  in  satisfying  the  objectives  listed  in  Table  ..1, 
and  to  determine  the  consistency  of  orbital  measurements  with  the  geostation¬ 
ary  data  taken  in  Prairie  Smoke  1  (and  therefore  determine  if  changes  should 
be  made  in  the  configuration  used  for  future  geostationary  experiments). 

2.2  Objectives  of  Prairie  Smoke  II  (U) 

(S)  Prairie  Smoke  II  represented  the  first  opportunity  to  mount  both  experi¬ 
ments  (geostationary  and  orbital)  simultaneously.  The  objectives  of  the 
Pr.'.irie  Smoke  II  test  series  were  primarily  to  obtain  sufficient  quality  and 
quantity  of  measurements  for  confirming  Prairie  Smoke  I  and  lb  results,  to 
provide  additional  information  concerning  the  horizontal  extent  and  altitude 
distribution  of  the  ASF  irregularities,  and  to  demonstrate  the  fading  depen¬ 
dence  on  frequency  by  taking  data  at  400  MHz  during  the  orbital  experiment. 

The  major  improvement  in  data  quality  planned  for  Prairie  Smoke  II  includes 
the  capability  of  real-time  digital  data-logging  for  the  geostationary  experi¬ 
ment  and  the  capability  of  A/D  conversion  of  the  orbital  data,  as  well  as 
data  collected  during  Prairie  Smoke  I  and  lb. 
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3.  EXPERIMENTAL  CONFIGURATIONS  (0) 


,s)  iiuring  Prai  rie  Smoke  .b  on  February  14-17,  1972.  observations  were  made 
of  the  scintillations  of  signals  radiated  fro.  US  Navy  navigations,  satel¬ 
lites  transmitted  at  about  ISO  Milt  and  moving  in  a  nearly  polar  orbit.  Since 
the  periods  of  observation  of  these  satellites  was  limited  t> 
considerations,  it  was  not  possible  to  use  a  single  receiving  site  for  all 
four  days.  It  was  then  assumed  that  the  greatest  intensity  of  fluctuation-. 
uould  be  found  when  the  satellite  bOS  passes  through  a  point  at  300-km  alti¬ 
tude  directly  over  Platteville.  For  the  passes  nearest  to  overhead,  the 
sate,  lire  LOS  was  traced  through  the  300-km  point  to  the  surface  of  the  earth. 
During  the  satellite  pass,  the  intersection  with  the  earth  weeps 
N-S  direction.  The  points  for  each  day  at  which  it  intersects  Interstate 
Route  80  were  determined,  and  established  as  the  observation  point  lor  each 
day  A  map  showing  these  site  locations  with  respect  to  the  transmitter  near 
Platteville  is  provided  in  Figure  3.1.  Subsequent  analysis  of  the  satellite 
passes  for  each  location  (Table  3.1)  showed  that  the  LOS  passes  over  Platte- 
vi He  at  altitudes  of  284  to  316  km,  in  good  agreement  with  the  predretions. 
The  availability  of  spaced  antenna  data  for  each  pass  is  provided  in  Table  3.2. 
(S)  Since  the  highly  mobile  receiving  station  for  the  orbital  expenme 
(subsequently  used  for  Prairie  Smoke  II)  was  not  yet  available  at  the  time-of 
Prairie  Smoke  lb,  a  station  was  equipped  by  rack -mounting  the  receivers  and 
chart  recorder  used  in  the  geostationary  experiment  of  Prairie  Smoke  I  m  a 
utility  van.  Additional  equipment  included  a  Sony  Type  366-4  4  cha 
deck  for  recording  the  product  detector  BFO  output  of  the  receivers.  Dual 
circularly-polarized  turnstile  antenn^nd  battery-operated  Vanguard  dual- 
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LOS  over  transmitter  for  orbital  passes  during  Prairie  Smoke  lb  (U) 
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Availability  of  data  for  orbital  masses  during  Prairie  Smoke  lb  (U) 
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gate  MOSF'HT  RI-  converters  were  used  for  receiving  the  signals  from  the  or¬ 
biting  satellites.  Data  was  subsequently  analyzed  by  playing  the  signals 
back  through  a  Hewlett-Packard  310A  wave  analyzer,  acting  as  a  tracking  fil 
ter  with  a  200 -Hz  bandwidth,  and  recorded  on  a  thermal-writing  recorder. 

(S)  During  Prairie  Smoke  II  on  April  24  -  May  12,  1972,  both  a  geos c at  ionary 
and  orbital  experiment  were  conducted.  The  geostationary  mobile  statior  was 
deplo>ed  near  Lance  Creek,  Wyo.  ,  the  more  southerly  site  location  used  during 
Prairie  Smoke  1.  Since  this  facility  required  at  least  one  full  day  for 
changing  site  locations,  a  single  site  was  used  for  each  day  even  though  the 
height  of  reflection  of  the  heater  changed  during  the  experiment. 

(^)  The  same  equipment  was  used  for  the  geostationary  experiment  as  was  used 
during  Prairie  Smoke  I.  Additionally,  the  product  detector  BFO  output  of  the 
receivers,  including  one  monitoring  WWV,  was  continuously  recorded  with  a  4- 
channcl  tape  deck.  Simultaneously,  the  receiver  detector  output  was  recorded 
on  a  Metrodata  model  616  digital  data-logging  system  sampling  at  2-^/2  samples 
per  second  on  each  of  three  channels.  Use  of  a  Metrodata  622  S  digital  tape 
reader  interfaced  with  a  PDP-Se  digital  computer  offers  the  potential  of  com¬ 
pact  storage  and  prompt  analysis  of  large  quantities  of  data.  Development  of 
the  software  to  implement  this  digital  data-vecovery  system  is  nearly  com¬ 
plete.  Real-time  records  were  also  obtained  by  recording  the  detector  output 
on  a  dual-channel  chart  recorder.  The  availability  of  the  geostationary  re¬ 
cords  is  given  in  Table  3.3. 

(S)  The  mobile  orbital  satellite  receiving  station  for  Prairie  Smoke  II  con 
sisted  of  a  23-ft  Winnebago  vehicle  shell  with  a  5-kW  gasoline  generator  as 
a  power  supply,  providing  a  mobile  facility  capable  of  traveling  between 
widely  separated  site  locations.  The  data  was  recorded  on  magnetic  tape  in 
the  manner  described  for  Prairie  Smoke  lb.  An  A/P  converter  with  an  8-channel 
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Table  3.3 


Availability  of  geostationary  data  during  Prairie 

GEOSTATIONARY  DATA  AVAILABLE 


DATE  (1972) 

TIMES  OF  RECORD  (to ST ) 

24  APRIL 

1247-1705 

25  APRIL 

0559  -  1200 

04  MAY 

0859  -  1945, 
2015-2206 

05  MAY 

1055-1701 

08  MAY 

0800-  1 143  *  * 

09  MAY 

0554-1118, 

1130-1501 

10  MAY 

1620-2202 

II  MAY 

1552-2200 

12  MAY 

1 102  - 151 1 

*  N-S  SPACING  ALWAYS  305  m 
**  TRANSMITTER  OFF 
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90 

90 
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90 

90 

122 

122 
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multiplexer  is  presently  available  for  use  with  the  PDP-8e  computer  to  per¬ 
form  more  rapid  data  reduction  from  the  tape  recordings  using  Fourier  and 
correlation  techniques  than  is  possible  by  visual  reduction  of  chart  recordings. 
(S)  Site  locations  for  the  orbital  experiment  were  chosen  so  that  the  LOS 
passed  near  the  300-km  point  over  the  transmitter  within  the  constraints  of 
time  to  change  site  locations,  transmitter  schedules  for  the  Platte”ille 
facility,  and  difficulty  in  crossing  the  continental  divide  in  the  months  of 
\pril  and  May.  It  was  also  desired  to  record  as  many  passes  as  possible  (even 
the  more  oblique  passes),  in  order  to  place  limits  on  the  size  of  the  volume 
comprising  the  modified  region  by  detecting  the  reduced  effect  at  the  edge 
of  the  region.  The  site  locations  for  the  passes  recorded  during  Prairie 
Smoke  II  are  given  in  Table  3.4.  The  distance  of  closest  approach  cf  the 
LOS  of  these  passes  from  the  300-km  point  over  Platteville  is  presented  in 
Table  3.5.  A  map  showing  the  location  of  these  sites  with  respect  to  the 
transmitter  near  Platteville  is  provided  in  Figure  3.2.  The  availability  ol 
spaced-antenna  data  for  the  orbital  passes  is  presented  in  lable  5.6. 
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Table  3.4 


Site 

DATE  (1972) 

locations  Jor  orbital 

TIME  (MST) 

experiment  during  Prairie  Smoke  11  (in 

SITE  LOCATION  AND  COORDINATES 

24  APRIL 

1336 

ESTES  PARK,  COLO.  (40.38°N,  I05.52°W) 

24  APRIL 

1409 

25  AMIL 

0856 

04  MAY 

0930 

BRUSH,  COLO.  (40.19°  N,  103.38°  W  ) 

04  MAY 

1606 

05  MAY 

1336 

WIGGENS,  COLO  (  40.09°  N  ,  104. 07°  W) 

08  MAY 

1300 

09  MAY 

0832  * 

HILLSDALE,  WYO.  (  4I.I7°N  ,  104.47°  W  ) 

09  MAY 

0900 

10  MAY 

2032 

KIMBALL,  NEB.  (4I.23°N,  103.72°  W  ) 

10  MAY 

21 10 

II  MAY 

2021 

ESTES  PARK,  COLO.  (40.38CN,  I05.52°W) 

12  MAY 

1257 

12  MAY 

1454* 

*TRANSM  TTER  OFF 
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Table  3.6 

Availability  of  orbital  data  during  Prairie  Smoke  II  (U) 


DATE  (1972) 

OBSERVATION 

LOCATION 

NUMBER  OF 
150  MHz 

ANTENNAS 
400  MHz 

24  APRIL 

ESTES  PARK,  COLO. 

4 

1 

25  APRIL 

4 

1 

04  MAY 

BRUSH,  COLO. 

4 

1 

05  MAY 

WIGGENS,  COLO. 

4 

1 

08  MAY* 

3 

0 

09  MAY 

HILLSDALE,  WYO. 

3 

1 

10  MAY 

KIMBALL,  NEB. 

2 

2 

II  MAY 

ESTES  PARK,  COLO. 

3 

2 

12  MAY 

3 

2 

‘TRANSMITTER  OFF 
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4.  RESULTS  0I;  PRAIRIE  SMOKE  lb  (U) 


(S)  The  satellites  of  the  IJS  Navy  Satellite  Navigation  System  arc  all  in 
approximately  circular  polar  orbits  at  heights  of  about  1025  km  above  the 
earth's  surface,  Since  their  times  of  passage  are  distributed  (though  some¬ 
what  irregularly)  throughout  the  day,  they  are  particularly  suitable  for 
transmission  studies  of  specific  regions  of  the  ionosphere  such  as  that  over 
Pl.U  tevi 1 le.  Unfortunately,  only  one  or  two  a  day  pass  sufficiently  nearly 
overhead  1’lattevillc  to  permit  their  use  as  a  consistent  diagnostic  for  a 
given  site  location,  and  these  may  not  be  at  suitable  times  on  a  given  day. 
However,  between  February  14-17,  1972  the  four  passes  described  in  Tables  5.1 
and  5.2  were  recorded. 

(S)  Figures  4.1  and  4.2  show  spaced  receiver  signals  recorded  during  a  por¬ 
tion  of  two  of  these  passes.  Subsequent  to  Prairie  Smoke  lb  it  has  become 
possible  to  play  the  tape  recordings  back  through  the  wave  analyzer  and  an 
A/D  converter  interfaced  with  a  PDP-8e  computer  as  described  in  Section  2  to 
provide  values  of  the  scintillation  index  over  4-second  intervals.  The  ex 
pression  for  the  scintillation  index  (S)  is  given  by 


2  1^2  1^  2  1^  2 

S  =  (ij  A  -  (  ±  Z  A  )  ),  (  ijA) 
k  N  j  n  N  i  n  N  j  n' 


(4.1) 


where  is  the  received  signal  amplitude  for  the  sample  and  N  is  the 
total  numbe1"  of  samples  taken  over  the  given  time  interval.  The  results  of 
that  analysis  is  presented  for  the  Kimball  and  Buford  pa' ses  in  Figures  4.^ 
and  4.4.  These  two  passes  are  analyzed  in  detail  here  since  they  represent 
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KIMBALL, NEB.  15  FEB.  72 


VVA^\/W~A^AvA 

SOUTH  ,  _ I - 


NORTH 


2149:45 


2149:50 


TIME  (UT) 


Figure  4.1 

Chart  records  showing  sample  of  orbital  data  (U) 
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BUFORD,  WYO.  16  FEB.  72 


Figure  4.2 

Chart  records  showing  sample  of  orbital  dat;  U) 
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Fading  intensity 
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BUFORD,  WYO 

4  sec  INTERVALS 


_ I _ 

1901*00 


TIME  (MST)  16  FEB  72 

Figure  4.4 

Fading  intensity  for  an  orbital  pass 
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the  two  passes  when  the  deepest  fading  was  observed  and  the  LOS  was  nearest 
t  ing  parallel  to  the  magnetic  field.  The  scintillation  indices  for  each 
ol  the  four  passes  when  the  deepest  fading  was  observed  are  presented  in 
iahle  4.1  with  an  estimate  of  the  angle  between  the  LOS  and  magnetic  field. 
Ihe  chart  recordings  at  these  times  are  shown  in  Figure  4.5.  The  transmitter 
mode  during  each  pass  is  provided  in  Table  4.2.  It  is  planned  to  perform  a 
complete  correlation  analysis  on  these  data  to  provide  size,  orientation,  and 
height  distribution  of  the  ASF  irregularities.  However,  for  the  antenna  spa- 
ongs  used  (76  m  at  Buford  and  122  m  at  Kimball)  along  a  N-S  baseline  the 
time  shift  for  maximum  cross  correlation  can  be  measured  by  comparing  chart 
recordings  run  at  faster  speeds  than  those  presented  in  the  following  figures, 
fhis  procedure  was  used  for  determining  the  height  of  the  irregularities, 
which  for  a  plane  earth  is  given  by 


V.  t 

2i  =  R  sin  (e)/(l  -  — — ) 


(4.2) 


where 

K  =  Range  to  satellite  from  observer 
c  =  Elevation  of  satellite  from  observer 
Vs  =  Satellite  velocity  parallel  to  baseline 
t  =  Time  shift  on  N-S  baseline 
d  =  N-S  antenna  spacing 

z.  =  Height  of  irregularities  above  tangent  plane. 
zi  converted  to  heights  above  the  earth's  surface  by 

h  =  (zr  +  zp)/  cos  yc  -  zr  (4.3) 
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SINCLAIR,  WYO. 
14  FEB.  72 


17  FEB.  72 

.  I  sec 

i - 1 


_ L _ 

0045:50 

TIME  (UT) 


Figure  4.5 

Sok”  ?b“;f  Sh0Win8  maximun'  fadlng  °bse™‘1  <^i"8  Prairie 
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who  re 

h  =  the  geocentric  height, 

2  =  the  radius  of  the  earth 

r 

y  the  angle  at  the  earth's  center  subtended  by  the  satellite  and  the 
point  of  observation  and  where  it  has  been  assumed  that  the  earth 
is  spherical. 

(S)  IVhen  the  records  for  the  Kimball  and  Biiford  passes  were  examined  and  the 
time  shifts  for  individual  fades  measured,  it  was  found  that  the  heights  cal 
ciliated  were  distributed  through  a  height  range  of  approximately  100  km. 

For  both  passes  a  height  of  225  km  was  most  frequently  measured  with  heights 
extending  up  to  the  F  layer  and  down  to  somewhat  less  than  200  km.  ihese 
height  distributions  are  shown  in  Figures  4.6  and  4.7. 

(S)  For  the  Kimball  pass  the  heights  observed  were  subdivided  into  those 
observed  in  two-second  intervals.  Figure  4.8  shows  these  height  ranges  on 
the  LOS  froin  the  observation  site.  The  mean  heights  observed  for  each  oi 
these  intervals  are  located  nearly  directly  over  the  transmitter  for  the  deepe 
fading  observed,  that  part  of  the  pass  shown  in  Figure  4.5. 

(S)  The  average  spatial  quasi -period,  a,  has  been  calculated  using 


a  =  (1  -  z-^R  sin  e)Vst 


(4.4) 


where  z. ,  R,  e,  and  V  have  been  previously  defined  and  t  is  the  lade  period 
at  the  observation  site.  Figure  4.9  shows  the  frequency  distribution  of  the 
quasi-periods  measured  at  Kimball.  For  a  Gaussian  auto-correlation  function 
(p)  of  the  form 


p(x)  =  e 


2  .  ,2 

-x  /d 


(4.5) 


SECRET 


SECRET 


(— 

co 

2 

O 

o 

o 

m 

5 

i 

co 

<n 

^3" 


< 

CD 

2 

ic: 

d 

UJ 

> 

oc 

UJ 

to 

03 

O 

CO 

I— 


X 

CO 

UJ 

2 


co 

1 

2 


O 

O 

Q 

UJ 

O 

< 

0. 

CO 

co 

< 

2 

2 

UJ 

I— 

2 

< 


Ul 

O 

>- 

o 

2 

UJ 

3 

O 

UJ 

oc 


10 


rt 

U 


.4 


Figure  4.6 

F-'rcquencv  distribution  of  h  for 


QUASI -PERIOD  OBSERVED 
KIMBALL,  1449:16-1450:00  MST 


Figure 

Frequency  distribution  of  a 
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the  structure  size  (d)  is  given  by 

d  =  a/2.56  (1  <>) 

['he  data  analysis  presently  planned  will  determine  the  validity  of  (5)  in 
describing  the  observed  fading. 

(S)  It.  is  interesting  to  compare  the  orbital  quas  i -peri  ods  of  i'igur  4.9  with 
the  corresponding  structure  sizes  found  in  Prairie  Smoke  1  by  the  geostation¬ 
ary  technique  (Bowhill  and  Mendenhall,  1972,  pp.  52-55).  The  latter  experiment 
showed  spatial  quas i -periods  between  200  and  600  m,  with  an  average  of  about 
500  m,  in  good  agreement  with  the  orbital  data  of  Figure  4.9. 

(S)  At  Buford  three  antennas  spaced  in  a  right  triangle  were  recorded  (N-S 
spacing  of  250  ft  and  N,  W-E  spacing  of  255  ft).  Since  the  time  correlation 
observed  at  the  three  antennas  was  almost  perfect,  the  pattern  on  the  ground 
appears  to  be  highly  elongated.  The  locus  of  points  of  maximum  cross  cor¬ 
relation  form  a  straight  line  known  as  the  line  of  maxima,  and  the  orientation 
of  the  line  of  maxima  is  entirely  dependent  on  the  satellite  look  angles  as 
demonstrated  earlier  by  Bowhill  and  Mendenhall  (1972). 

(S)  The  geometry  for  the  Buford  site  location  and  a  pass  which  is  east  of 
overhead  is  given  by  the  following  relations: 

9  =  if  +  arctan  (sin  5/(tan  6  cot  e  -  cos  6))  (4.7) 

where 

if  =  (180  -  8)  where  8  is  the  satellite  azimuth  from  the  observation  site 
6  =  (D  +  if)  where  D  is  the  magnetic  dip  in  the  region  of  the  irregularities 
G  =  Dip  of  the  magnetic  field  in  the  region  of  the  irregularities 
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t  =  Elevation  of  the  satellite  from  the  observation  site. 

1-or  the  antenna  configuration  used  at  Buford  the  observed  orientation  is  given 
by  the  following  relations: 

0  =  arctan  (d^t/dtj)  (1.8) 

who  re 

dj  =  N,  h-1;  spacing 

tj  Time  shift  for  maximum  cross  correlation  between  the  N,  W  and  I  an¬ 
tennas 

d  =  N,  W-S  spacing 

t  =  Time  shift  for  maximum  cross  correlation  between  the  N,  IV  and  S  an¬ 
tennas  . 

A  few  samples  of  data  were  analyzed  by  measuring  the  time  shifts  for  maximum 
cross  correlation  from  chart  recordings.  The  results  are  shown  in  Table  4.3, 
where  the  magnetic  dip  was  taken  to  be  14  deg  for  the  calculated  values  based 
on  the  geometry  of  the  satellite  pass. 

(S)  The  agreement  between  the  measured  orientation  in  Table  4.3  and  the 
orientation  calculated  on  the  basis  of  field-alignment  of  the  irregularities 
confirms  that  the  assumption  of  field-alignment  is  well  justified. 
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S.  RESULTS  OF  PRAIRIE  SMOKE  II  (IJ) 

(S)  m  Prairie  Smoke  II  for  the  first  time,  orbital  and  geostationary  data 
were  available  simultaneously. 

5.1  Results  of  the  Geostationary  Experiment  (U) 

(S)  Based  on  the  geostationary  data  listed  in  Table  3.3,  the  scintillation 
index  S  has  been  calculated  using  the  formula  (1)  of  Section  3.3.  It  was 
found  to  vary  greatly  during  the  experiment.  To  facilitate  comparison  with 
other  data,  it  has  been  calculated  for  a  2-minute  interval  fifteen  minutes 
.liter  each  hour  when  the  transmitter  was  operating.  These  results  are  given 
in  lable  5.1.  Comparison  of  these  data  with  parameters  such  as  the  data  from 
the  phased  array,  the  spread  observed  on  the  lonosonde,  the  heater  reflection 
.height  and  the  heater  frequency  relative  to  the  critical,  will  give  insight 
into  a  yield  model  of  the  disturbed  region. 

(S)  An  example  of  the  deep  fading  frequently  observed  is  shown  in  Figure  5.1. 
In  this  case,  the  autocorrelation  function  was  found  to  have  an  e-folding 
time  of  about  11  seconds  (Figure  5.2).  Its  shape  generally  supports  the 
Gaussian  form  assumed  for  the  autocorrelation  function  in  Equation  4.6. 

(S)  A  survey  of  the  real-time  chart  recordings  indicates  that  the  fading  was 
generally  deeper  after  1800  MST.  Typical  values  of  S  during  the  morning  and 
afternoon  hours  were  1-5%,  while  deep  fading  (So,  5-25%)  was  frequently  observed 
during  the  evening  hours.  The  results  of  the  orbital  experiments  have  not 
shown  the  same  dependence  on  time  of  day,  suggesting  that  at  least  part  of  the 
day /night  variation  is  an  effect  of  the  change  in  height  of  the  ionosphere 
from  day  to  night.  Apparently,  the  disturbed  region  was  closer  to  the  satel- 
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Table  5.1 

Scintillation  Index 


Cal culated 

From  Geostationary  Satellite 

Data  (U) 

DATE 

TIME  (MST) 

SCINTILLATION 

19  72 

INDEX  (%) 

24  April 

11 :  15 

<  1.0 

12:  15 

<  1.0 

13:15 

4.3 

14:15 

2.7 

15:15 

2.2 

16:15 

<  1.0 

25  April 

6:15 

<  1.0 

7:15 

'  3.3 

8.  15 

4.5 

9:15 

5.6 

10:15 

<  1.0 

11:15 

<  1.0 

4  May 

9:15 

3.3 

10:15 

2.2 

11:15 

3.6 

12:15 

<  1.0 

13: 15 

1.7 

14:15 

2.5 

15: 15 

<  1.0 

16:15 

2.8 

17: 15 

4.8 

18:15 

9.6 

19:15 

5.7 

20:15 

12.0 

21:15 

10.0 
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DATE 


5  May 


9  May 


10  May 


1 1  May 


1 2  May 


Table  5.1  cont. 


SCINTILLATION 
_ INDEX  (%) 


11 :  15 
12:15 
13:15 
14:15 
15:  15 
16:15 
6:  15 
7:  15 
8:  15 
9: 15 
10: 15 
11:15 
12:15 
13: 15 
14:15 
17:  15 
18:  15 
19:15 
20:15 
21: 15 
16:15 
17:15 
18:15 
19:15 
20: 15 
21:15 
11:15 
12:15 
13:15 
14:15 


<  1.0 
3.1 
1.0 

3.3 
2.9 
7.5 

2.4 
1.2 
3.2 

2.5 
4.0 

2.6 

3.1 
3.9 

<  1.0 
<  1.0 

4.1 

4.4 
5.3 

25:5 

3.5 

2.9 

4.9 

<  1.0 

5.7 
13.7 
<  1.0 
<  1.0 
<  1.0 
<  1.0 
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LANCE  CREEK ,  WYO.  4  MAY  72 
SCINTILLATION  INDEX  =  16.1  % 


L 


2020 


2021 


2022 

TIME(MST) 


2023 


2024 


Figure  5 . 1 

A  sample  of  the  geostationary  data  (l!) 
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AUTOCORRELATION  COEFFICIENT 
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Figure  5.2 

Autocorrelation  coefficient  for  a  sample  of  the  geostationary  data.  (U) 
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lire  LOS  during  the  ni^ht  hours.  In  fact,  orbital  data  seemed  to  show  little 
dependence  of  yield  on  time  of  day.  The  maximum  value  of  S  for  a  4-second 
interval  was  found  to  be  12.9%  at  Pine  Bluffs,  Wyo.  on  Dec.  6,  1971  and  12.6% 
at  Kimball,  Neb.  on  Feb.  15,  1972.  These  were  two  afternoon  passes.  In  Sec¬ 
tion  5.2  an  evening  pass  on  May  10,  1972  is  described  showing  a  similar  value 
of  13.3%. 

(S)  The  frequency  distribution  of  the  spatial  quasi-period  observed  at  the 
geostationary  site  is  shown  in  Figure  5.3  for  a  sample  of  the  data  collected 
on  May  10,  1972.  It  shows  general  agreement  with  the  results  of  the  orbital 
experiment  and  with  the  previous  geostationary  data  from  Prairie  Smoke  1. 

(SJ  A  survey  of  the  real-time  cnart  records  was  made  to  determine  onset  and 
decay  times  and  any  dependence  with  transmitter  power  and  fof'2/f^.  Hie  scat¬ 
ter  plots  in  Figure  5.4  show  the  results  for  the  onset  times.  The  decay  times 
were  observed  to  vary  from  a  few  minutes  to  in  excess  of  a  half-hour  and  were 
gererally  longer  at  night.  No  very  pronounced  correlation  is  found  with 
heater  frequency  or  input  power. 

(S)  Fxamination  of  the  geostationary  records  show  that  the  drift  velocity 
was  typically  30-60  m/sec  and  moving  toward  either  the  east  or  west.  On  May  11, 
1972  at  2000  MST  the  direction  of  the  structure  motion  was  observed  to  grad¬ 
ually  shift  from  being  easterly  to  westward,  as  found  in  Prairie  Smoke  1. 

(S)  The  geostationary  data  is  being  further  analyzed  to  provide  information 
concerning  the  axial  ratio,  orientation,  and  drift  velocity  based  on  the  cross - 
correlation  functions  of  spaced-antenna  records. 

5.2  Results  of  the  Orbital  Experiment  (U) 

(S)  Only  a  limited  amount  of  the  orbital  data  has  been  fully  analyzed  at  this 
time,  but  reported  in  this  section  are  the  results  of  a  pass  observed  from 
Kimball,  Neb.  with  spaced-antenna  records  were  of  good  quality  for  both  the 
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Figure  5.3 

Spatial  quasi -periods  observed  at  the  geostationary  site.  (U) 
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Figure  5.4 

Onset  times  observed  at  the  cationary  site.  (U) 
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150  and  400  MHz  signals.  The  index  S  has  been  calculated  at  4-second  inter¬ 
vals  and  is  shown  in  Figure  5.5.  A  portion  of  the  spaced-antenna  chart 
recording  for  the  400 -MHz  signal  is  shown  in  Figure  5.6. 

(S)  When  the  data  of  Figure  5.5  are  used  to  form  a  scatter  plot  (Figure  5.7) 
it  is  evident  that  the  ratio  of  the  fading  depth  (while  showing  considerable 
variability)  corresponds  approximately  to  an  inverse  frequency  dependence  for 
the  scintillation  index.  This  result  agrees  with  the  diffraction  theory  of 
Briggs  and  Parkin  (1963)  for  irregularities  located  in  the  far-field  or  Fraun¬ 
hofer  region.  The  conditions  for  this  result  are  that  the  RMS  phase  deviation 
be  less  than  unity  and  that 


Az/D  >  tt/3 


(5.1) 


where 

z  is  the  effective  distance  from  the  scatterers,  given  by  ZjZ.,/(Zj  +  zj 
Zj,  z2  are  the  distances  from  the  irregularities  to  the  observer  ar.d 
source,  respectively 

A  is  the  wavelength  of  the  source  signal 

IJ  is  the  e-folding  distance  for  the  transverse  correlation  function  of 
the  irregularities. 

(S)  The  height  distribution  of  the  irregularities  has  been  calculated  in  a 

I 

manner  similar  to  that  described  in  Section  4  and  is  shown  in  Figure  5.8  for 
both  frequencies.  The  most  striking  feature  of  this  variation  of  height  with 
time  is  the  rapid  increase  at  2032:10  MST  and  the  corresponding  abrupt  de¬ 
crease  at  2052:55.  It  is  easy  to  explain  this  variation  if  one  considers  the 
geometry  of  the  satellite  pass. 
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KiMBALL  ,  NEB.  10  MAY  72 
TIME  CONSTANT  =  5  msec 
FREQUENCY  =  400  MHz 
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Figure  5.6 

Sample  of  spaced  antenna  records  for  an  orbital  pass  (U) 
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Figure  5.7 

Scintillation  ind^x  observed  at  two  frequencies  (U) 
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(S)  Imagine  that  the  disturbed  region  is  confined  within  a  cylinder  parallel 
to  the  lines  of  the  earth's  magnetic  field,  with  its  projection  on  a  300-km 
altitude  plane  a  circle  witn  radius  X.  Let  the  center  of  this  circle  be  a 
distance  to  the  magnetic  north  of  Platteville.  Then,  for  various  assump¬ 
tions  about  X  and  xQ ,  heights  can  be  calculated  for  the  entry  and  emergence 
of  the  satellite  LOS  from  the  cylinder  containing  the  disturbed  region,  figure 
5.9  shows  these  curves  for  three  cases  of  interest.  The  area  contained  by 
these  curves  represents  the  height  range  over  which  scintillations  may  be  ex 
pected,  as  a  function  of  time  during  the  satellite  pass. 

(S)  Of  course,  a  particular  height  will  only  give  scintillations  as  long  as 
substantial  ionization  is  present.  Assume,  therefore,  that  200  km  is  the 
lower  limit  for  the  disturbed  region.  It  is  evident  that  the  model  for 
x()  =  300  km,  X  =  50  km  gives  a  height  for  the  scintillation  region  that  increases 
from  200  km  at  2032:00  to  200-450  km  at  2032:40;  decreasing  again  to  200  at 
2033:10.  This  is  exactly  the  height  variation  shown  in  Figure  5.8.  Further, 
since  the  depth  of  the  scintillation  increases  as  the  thickness  of  the  dis¬ 
turbed  region  increases,  the  scintillations  should  have  started  around  2052:00, 
passed  through  a  maximum  at  2032:40  and  disappeared  at  2033:10.  Comparison 
with  Figure  5.5  confirms  this  model.  The  additional  scintillation  seen  be¬ 
tween  2033:10  and  2033:30  must  arise,  according  to  Figure  5.8,  from  irregu¬ 
larities  below  200  km  altitude  at  the  southern  edge  of  the  region. 

(S)  The  quite  different  behavior  of  the  height  shown  in  Figure  4.8  is  easily 
explained  in  terms  of  the  different  geometry  involved.  Figure  5.10  is  a 
plot  similar  to  Figure  5.9  for  the  satellite  pass  of  Figure  4.C  and  Figure 
3.  Comparison  with  Figure  4.3  suggests  that  in  this  case  x^  =  -30  km, 

=  50  km;  and  because  of  the  different  shape  of  the  contour  on  figure  5.10, 
e  effective  height  increases  continuously  through  the  satellite  pass. 
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Figure  5.9 

Height  ranges  calculated  for  an  orbital  pass  (U) 
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Height  ranges  calculated  for 
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(S)  The  RMS  phase  fluctuation  A  for  a  wave  propagated  through  a  thickness 
Ah  of  a  random  medium  has  been  given  by  Briggs  and  Parkin  (1963)  as 


2  2 

7  TT_I)  E 

Y“  =  r-2 


1/2 


2  2  2  ; 
0  sin  ip  +  D  cos  \p 


r^  A^  AN^  Ah  sec  i 
e 


(5.2) 


where  H  and  I!  are  the  e-folding  correlation  distances  perpendicular  and  par¬ 
allel  t>  :  magnetic  field,  respectively;  ip  is  the  angle  between  the  wave 

7 

normal  and  the  magnetic  field;  re  is  the  classical  electron  radius;  AN“  is 
the  mean  square  electron-density  fluctuation;  and  i  is  the  angle  between  the 
wave  normal  and  the  vertical. 

(S)  Since  T  >>  I)  it  may  be  rewritten  in  more  familiar  notation  as 


2 

Y  = 


5/2 


f  4  7 

N,  ,AN 


Cj-)  (j^j— )  Ah  D  sec  i  cosec  ip 


(5.5) 


where  f'N  is  the  plasma  frequency  and  f  the  wave  frequency.  Since  f^  varies 
with  altitude,  this  may  be  rewritten  in  the  following  way: 


2 

Y 


5/2 


c-f 


D  sec  i  cosec  \p  f  (N/N  )  dh 

nr 


(5.4) 


on  the  assumption  that  the  fractional  RMS  electron-density  fluctuation  is  in¬ 
dependent  of  altitude.  The  integral  is  typically  about  120  km. 

(S)  Tor  the  scale  of  irregularities  corresponding  to  D  =  150  m,  far-field 
conditions  prevail,  and  the  RMS  amplitude  fluctuation  e  is  therefore  given  by 


e  =  y/*^2  . 


(5.5) 
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Figure  5.5  illustrates  an  orbital  pass  and  shows  that  e  goes  inversely  as 
frequency,  confirming  the  far-field  condition.  The  following  values  of  the 
parameters  apply  to  this  pass: 
y  =  0.075 
-  =  0.052 
f  =  400  MHz 
foF2  =  6.0  MHz 
i  =  26.6  deg 
i -  13.1  deg 
D  =  120  m 

/(N/N  )2dh  =  120  km 
m 

Solving  equation  (5.5)  with  these  values,  the  RMS  fractional  fluctuation  in 
electron  density  is  found  as  6.9  x  10  3. 
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(>.  DISCUSSION  OF  RESULTS  OF  PRAIRIE  SMOKE  lb  and  II  (U) 

(S)  The  orbital  data  from  Prairie  Smoke  lb  and  II  generally  confirmed  the 
conclusions  from  the  geostationary  data  of  Prairie  Smoke  I  and  the  orbital 
data  of  Prairie  Smoke  la.  Simultaneous  observations  at  150  and  400  Mile 
showed  an  inverse  frequency  variation  of  scintillation  depth,  and  confirmed 
the  Briggs-Parkin  diffraction  theory.  The  structures  of  ASF  were  confirmed 
to  be  field-aligned.  The  observations  indicated  a  scattering  model  with 
striations  contained  in  a  field-aligned  cylinder,  radius  50  km,  with  its  lower 
face  at  200-km  altitude,  and  centered  approximately  over  the  heater  at  300-km 
altitude.  The  electron-density  fluctuation  was  about  0.1  percent  RMS  at  the 
center  of  the  region. 

(S)  Continuous  geostationary  scintillations  showed  marked  fluctuations  in 
scintillation  index  which  may  be  correlated  with  changes  in  height  of  the 
F2  layer. 

(S)  Based  on  these  data,  a  scattering  model  has  been  prepared  for  preliminary 
applications  evaluation  and  is  described  in  Sections  7  and  11. 
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7.  PURPOSE  AND  SCOPE  OP  MODEL  (U) 

(S)  In  order  to  evaluate  the  possible  applications  of  irregularities  produced 
by  heating  of  the  ionospheric  F  region,  it  is  necessary  to  describe  the  ir¬ 
regularities  by  a  simple  scattering  model.  In  determining  the  scope  of  sucli  a 

1  ^  1  ' -s  useful  to  consider  two  boundary  conditions  that  circumscribe  its 

degree  of  elaboration: 

(  l .)  It  should  not  be  simply  a  compilation  of  cross  sections  observed 
between  pairs  of  transmitting  and  receiving  stations  in  different 
locations  at  different  frequencies,  but  should  be  a  self-consistent 
description  of  the  scattering  medium  itself;  albeit  in  simplified 
terms . 

(iij  It  need  not  be  devised  to  be  in  agreement  with  any  particular  aspect 
of  plasma  instability  theory,  unless  that  agreement  has  been  fullv 
verified  empirically. 

(S)  Evidently,  this  approach  will  lead  to  an  empirical  description  of  the 
scattering  irregularities  which  could  be  aptly  described  as  an  "engineering 
model  rather  than  a  "physical  model".  Nevertheless ,  given  the  present  un¬ 
certainties  in  the  physical  mechanisms  involved,  it  seems  to  be  the  only  ap 
proach  that  can  lead  to  a  description  of  the  disturbed  region  with  a  relatively 
small  number  of  parameters. 

(S)  The  model  should  be  capable  of  estimating  numerically  the  following 
quantities: 

(a)  The  total  power  scattered  in  a  given  direction  by  the  disturbed  region 
when  illuminated  by  a  given  frequency  from  a  transmitting  antenna 
with  any  prescribed  polar  diagram  and  pointing  direction  relative  to 
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the  center  of  the  disturbed  region. 

(b)  The  distribution  in  range  of  the  scattered  power  (for  multipath 
considerations) . 

(c)  The  amplitude  and  phase  fluctuation  to  be  expected  in  a  signal  of 
specific  frequency  transmitted  through  the  irregular  region  from  the 
ground  to  space  (or  vice  versa). 

(S)  Because  of  the  variability  of  the  phenomenon,  it  seems  clear  that  some 
limitations  must  be  imposed  on  this  model  these  are  described  in  the  next 
sect  ion . 

7.1  Limit ations  of  the  Model  (Li) 

(S)  This  first  scattering  model  ^s  limited  to  the  following  conditions: 

(i)  Heating  from  the  Plattcvillc  transmitter,  with  full  power,  the 
antenna  pointed  directly  overhead. 

(ii)  Ordinary-mode  heating. 

(iii)  Heater  frequency  less  than  foF2. 

(iv)  Daytime  conditions  typical  of  the  larger  cross  sections  and  scintil¬ 
lations  observed. 

(v)  On-frequency  scatter  only  (rather  than  plasma  line). 

(S)  It  is  possible  that  optimum  conditions  may  later  be  discovered  giving 
much  greater  yields  than  are  indicated  by  the  present  model.  Nevertheless, 
it  gives  a  basis  on  which  to  build  a  later  model. 

(S)  The  model  is  treated  in  three  parts,  each  in  a  separate  section:  the 
magnetic  field;  the  on-frequenty  scatter  (OFS) ;  and  the  artificial  spread  F 

(ASF)  . 
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8.  MAGNETIC  El  ELI)  MODEL  (II) 


(S)  Both  OL'S  and  ASE  have  been  established  as  critically  dependent  on  the 
magnetic  field  direction.  Thorne  (1972)  has  used  a  99 -term  magnetic  field  model 
to  predict  the  properties  of  the  magnetic  field  over  Platteville.  A  simpler 
approach  is  suggested  here  which  apparently  gives  adequate  accuracy  for  most 
purposes . 

(U)  The  magnetic  dip  6  around  a  magnetic  dipole  is  given  by  the  equation 


tan 


2  tan  X  n 
mO 


(8.1) 


and  is  constant  with  altitude  over  a  point  P  with  magnetic 
Therefore,  the  actual  magnetic  dip  over  another  point  Q  in 
of  I1  and  on  the  same  magnetic  meridian  can  be  conveniently 
equation 


latitude,  A  .  . 

mO 

the  neighborhood 
represented  by  the 


tan  6=2  tan  X  .  (8.2) 

m 

The  magnetic  latitude  X  of  the  point  Q  is  given  to  a  close  approximation  by 
the  equation 


X  =  X  .,  +  (A  -  A-Jsec  a  (8.3) 

m  ml  0 

where  X and  X  are  the  geographic  latitude  of  P  and  Q,  respectively,  and  a  is 
the  magnetic  declination. 

SECRET 

54 


*vl 


SECRET 


(S)  Using  these  equations,  the  altitude  h  of  the  specular  reflection  point 
may  he  calculated  for  any  bistatic  radar  pair.  In  addition,  the  magnetic  dip 
6^  at  the  point  P  may  be  adjusted  so  as  to  fit  the  observed  altitudes  of  field 
aligned  scatter.  Figure  8.1  shows  a  comparison  of  h^  for  White  Sands  with  the 
Raytheon  calculation  of  Thome  (1972)  and  magnetic  dips  of  67.7  and  hS.2  deg 
over  Platteville.  The  effect  of  small  changes  in  magnetic  dip  by  far  outweighs 
the  small  differences  between  the  dip-latitude  model  and  the  99-coefficient 
model.  Accordingly,  the  dip  latitude  model  is  adopted  in  the  succeeding  scc- 
t ions. 
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9.  SMALL- IRREGULARITY  MODEL  (IJ) 

(S)  In  treating  on-frequency  scatter  (Ob'S),  there  are  several  questions  that 
need  to  he  resolved: 

(1)  How  ear  the  aspect  dependence  of  the  scattering  irregularities  be 
distingi  ished  from  the  effects  of  thickness  of  the  scattering  layer? 

( 2 )  How  can  the  scattering  cross  section  best  be  characterized  when  the 
scattering  medium  is  comparable  with  the  size  of  the  antenna  beam? 

(3)  To  what  extent  is  the  range  spread  of  the  irregularities  determined 
by  the  combination  of  aspect  sensitivity  and  disturbed-region 
thickness,  compared  with  the  effects  of  the  latitude  extent  of  the 
heated  region? 

These  three  questions  will  be  addressed  in  turn. 

y.l  Thickness  and  Aspect  (U) 

(S)  Imagine  that  the  scatter  irregularities  are  field-aligned,  m.J  are  dis¬ 
tributed  with  altitude  h  according  to  a  function  of  1(h)  given  by 

1(h)  =  exp (-h  -  h^) V A"-]  (y.l) 


where  h  is  the  altitude  of  reflection  of  the  heater,  defined  as  where  the 
r 

maximum  of  the  Airy  function  occurs;  and  A  is  the  e-folding  semi  thickness  of 
the  disturbed  region.  The  center  of  monostatic  scatter  from  the  irregularities 


will  be  at  the  height  of  h^  of  Section 


Let  the  scattering  polar  diagram  of 


tie  irregularities  in  the  elevation  plane  have  an  e-folding  semiwidth  6,  shown 


by  Booke’’  (1956)  to  be  given  by 
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0  =  A//2  tt  L  (9  ,  j  j 

Wh  re  is  the  wavelength  and  L  is  the  fcT  -folding  correlation  distance  along 
the  magnetic  field.  Scatter  from  a  height  h  rather  than  h^  will  have  a 
scattering  angle 


0  =  2(h  -  hq)/R  sin  5  (9..' 5) 

where  R  is  the  slant  range.  The  fraction  of  the  total  scattered  energy  which 
appears  in  the  angle  range  0  to  0  +  de  is  given  by 

J  (9) d9  -  (l//rT  0)exp(-02/0‘')d9  .  (9  .  >P 

Therefore,  the  fraction  J(h)dh  of  the  total  scatter  which  appears  to  come 
from  height  h  is  given  by 

J  (h)  dh  =  (2/vV  0R  sin  6)exp[-Ci  -  hQ)2/(R0  sin  6/2)  2]  .  (9.5) 

The  combined  effects  of  aspect  and  thickness  are  therefore  obtained,  at  a 
given  slant  range,  by  the  integral 

K  =  /I(h)J(h)dh  ,  (9.6) 

assuming  that  the  illuminating  antenna  gain  is  nearly  constant  over  the  height 
interval  involved. 
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(S)  Now  it  is  certain  that  the  effects  of  1(h)  play  a  dominant  part  in  de¬ 
termining  the  effective  cross  section;  so  that  term  cannot  be  neglected. 
Therefore,  until  further  experimental  evidence  is  available,  it  seems  most 
appropriate  to  make  the  assumption  that 

A  <<  RB  sin  6/2  .  (•’ •  ) 

Under  these  circumstances,  J (h)  reduces  to  a  6- function 

J  (h)dh  =  6  (h  -  h  )dh  (■-’ -81 


and 


K  =  I(h())  -  exp[-(hQ  -  hr)2/A2]  .  PS!n 

Justification  for  this  position  can  be  presented  as  follows:  Imagine  that  the 
irregularities  extend  all  the  way  through  the  disturbed  region,  so  that  L  -  A. 
What  would  the  value  of  A  need  to  be  so  that  he  inequality  (9 .  . )  became  an 
equality?  Inserting  Equation  (9.2)  into  (9.7)  the  equality  becomes 

A2  =  R>  sin  6/2/2tt  .  (  9‘ 10 J 

Using  R  =  1000  km  and  X  =  2m,  A  =  .5  km,  approximately.  It  therefore  follows 
that  if  the  irregularities  extend  through  the  entire  region,  and  are  longer 
than  1  km,  the  replacement  of  (9.5)  by  (9.8)  is  fully  justified.  The  Booker 
scattering  formula  then  gives  a  cross  section  per  unit  area,  the  height  de¬ 
pendence  having  been  integrated  out  of  (9.6)  to  give 
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K(h)  =  exp  [  -  (hq  -  hr)2/A']  .  (9-H) 

(S)  At  a  given  range,  therefore,  the  scattering  is  essentially  from  an  east- 
west  line  at  height  h  .  The  range  of  elevation  angles  of  the  radar  over  which 
the  scatter  is  observed  will  be  determined  by  the  convolution  of  the  antenna 
beam  with  (9.1);  the  latter  being  equivalent  to  a  Gaussian  distribution  with 
e-folding  semi angle  A/R  sin6.  Typically,  this  elevation  spread  is  determined 
by  the  antenna  polar  diagram,  so  the  only  information  available  at  a  given 
range  is  the  scattering  cross  section  and  the  altitude  h  . 

9 . 2  Cross  Section  Rationale  (U) 

(S)  As  indicated  in  Section  9.1,  scattering  at  a  given  range  comes  from  an 
arc  at  altitude  h^.  As  a  function  of  range,  therefore,  the  returns  must  come 
from  a  surface  determined  by  the  variation  with  range.  This  surface  is 
depicted  in  Figure  8.1,  and  its  altitude  depends  on  the  magnetic  dip.  There¬ 
fore,  it  is  meaningless  to  speak  of  a  cross  section  per  unit  volume,  in  that 
the  t'  ickness  in  scatter  of  volume  is  much  smaller  than  the  antenna  beam,  and 
probably  less  than  1  km.  Therefore,  it  is  best  to  speak  of  a  scattering  cros: 
section  per  unit  area  of  the  h^  surface.  Integrating  the  standard  Booker  ex¬ 
pression  with  altitude  at  a  given  slant  range  gives  the  following  expression 

for  cross  section  o  '  per  unit  area 

a 

oa  =  (r3ART2/AN4)exp(-8TT2T2/X2)  (9.12) 

where  AN/N  is  the  fractional  fluctuation  in  electron  density,  and  T  is  the 
/e-folding  correlation  distance  of  the  ionization  irregularities  transverse 
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to  the  magnetic  field.  The  model  proposed  here  suggests  that  the  variation 
of  electron-density  fluctuation  with  latitude  and  height  is  given  by 

CAN/.V)2  =  "c4N/Tt0  exp[-(X  -  X0)2/(X2/R02)  -  O'  -  lir>2/-n  O’-13) 


wheie  (AN/N)2  is  the  electron-density  fluctuation  across  the  magnetic  field 
at  the  height  of  reflection  hr  of  the  heater  free aency ,  at  the  latitude  XQ 
where  the  greatest  energy  is  deposited;  X  is  the  e-folding  distance  of  the 
heating  effect  in  the  north-south  direction,  and  RQ  is  the  radius  of  the  earth. 
The  exponential  factor  in  this  latter  equation  will  be  used  (with  h  =  hQ)  m 
the  next  subsection  to  illustrate  the  way  in  which  aspect  reacts  on  to  tins 

problem . 


9.5  Range  Spread  (U) 

(S)  The  question  arises  as  to  how  the  parameters  A,  X,  and  6Q  can  be  deduced 
from  the  experimental  data,  figures  9.1  through  9.4  illustrate  an  appropriate 
procedure.  They  show  the  relative  cross  section  as  a  function  of  slant  range 
for  a  heated  region  centered  30  km  north  of  Platteville,  the  magnetic  dip  of 
86.2  deg,  reflection  height  hr  of  275  and  300  km,  heating  radii  X  of  50  and 
100  km,  and  thicknesses  2A  of  10  through  200  km.  Comparison  of  Figures  9.1 
and  9.2  illustrate  how  the  range  extent  of  the  return  is  governed  at  small 
thicknesses  entirely  by  aspect,  but  at  larger  thicknesses  by  the  heated  region 
extent.  Figures  9.3  and  9.4  show  that  the  effect  is  even  more  striking  with 


X  =  50  km. 

(S)  This  information  may  be  presented  another  way,  shown  on  Figure  9.5.  For 
X  =  10u  km,  this  diagram  shows  the  importance  of  the  height  of  reflection  m 
determining  the  highest  apparent  cross  section  (as  a  function  of  range).  For 
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RELATIVE  CROSS  SECTION 
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SLANT  RANGE  (km) 


800  850  900  950  1000  1050 


LATITUDE  OF  REFLECTION  POINT  (deg) 


16.0  15.0  14.0  13.0  12.0 

ELEVATION  (deg) 

Figure  9.1 

Relative  cross  section  vs  slant  range  and  elevation  for  various  values  of 
2A  with  X  =  100  and  hr  ~275  (U) 
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SLANT  RANGE  (km) 
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LATITUDE  OF  REFLECTION  POINT  (deg) 
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Figure  9.2 

Relative  cross  section  vs_  slant  range  and  elevation  for  various  values  of 
2A  with  X  =  100  and  h  =  300  (U) 
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Figure  9.3 

Relative  cross  section  vs  slant  range  and  elevation  for  various  values  of 
2A  with  X  =  50  and  hr  =  275  (U) 
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RELATIVE  CROSS  SECTION 


SLANT  RANGE  (km) 


ELEVATION  (deg) 


Figure  9.4 

Relative  cross  section  vs_  slant  range  and  elevation  for  various  values 
2A  with  X  =  50  and  hr  =  300  (U) 


SE 


SECRET 


(6»p>  "j>  boj  aaniuvn 


(6»p)  UOJ  N0IJ.VA313 


Figure  9.5 

Changes  in  slant  range  and  elevation  for  maximum  cross  section  vs  h  for 
various  values  of  2A  (U) 
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example,  for  /A  =  20  km,  the  h  must  lie  between  260  and  a  10  kin  to  avoid  more 
than  20  dB  loss  relative  to  the  optimum  cross  section.  Varying  the  dip  ingle 
6  changes  the  magnetic  aspect,  as  does  changing  the  northern  offset  x()  oi  the 
most  intense  part  of  the  heated  region.  These  effects  have  been  shown  by  ad¬ 
ditional  horizontal  axes  for  hr  corresponding  to  two  values  of  5Q  and  xQ.  A 
similar  set  of  curves  is  shown  in  Figure  9.6  for  50  km  radius;  interestingly, 
only  the  thinner  regions  arc  affected. 

(S)  The  slant  ranges  of  the  greatest  cross  section  om  are  shown  on  Figures 
9.7  and  9.8.  It  is  interesting  to  note  that  the  thicker  the  region,  the  smal¬ 
ler  the  effect  of  heater  reflection  height  on  the  range  where  the  greatest 
cross  section  occurs.  This  is  probably  the  most  sensitive  way  to  determine 
phe  value  of  A. 

(S)  The  combined  effects  of  magnetic-field  geometry  and  heated-region  offset 
can  be  represented  by  the  following  equation  for  the  reflection  height  of  tin 
heating  wave  to  give  the  greatest  scattering  cross  section  at  White  Sands: 

h  =  ’285.6  -  17.1  (6Q-68 . 2)  -  .137  xQ  .  (9.14) 

(S)  The  thickness  of  the  heated  region  is  probably  sufficiently  small  that 
one  may  neglect  the  fact  that  its  thin  dimension  is  along  the  magnetic  Held 
axis  rather  than  in  the  vertical  direction.  This  is  not  the  case,  however, 
with  the  ASF  treated  in  Section  11. 

9.4  Additional  Work  Needed  on  this  Model  (U) 

(S)  The  following  procedure  is  suggested  as  a  means  of  arriving  at  the  ranges 
of  values  of  t..e  parameters  A,  X  and  x^.  The  reference  dip  <5q  can  be  dr 
termined  at  VHF  and  UHF  from  the  range-elevation  diagram  of  the  White  Sands 
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200  250  ~ 300  350  400  450 

hr(km),  A0  =  30km  N,  80=  67.7  deg 


200  250  300  350  400  450 

hr(km),  x0=30km  S,  80s  68.2  deg 

Figure  9.8 

Relative  value  of  the  maximum  cross  section  vs_  h  for  various  values  of 
XQ  and  6q  with  X  =  50  (U) 
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radar.  This  capability  is  not  available  at  frequencies  lower  than  1.10  Ml  lz , 
fortunately,  it  may  be  assumed  that  exactly  the  same  6  applies  at  those  fre¬ 
quencies.  The  range  of  <5^  needs  to  lie  established  as  it  is  important  to  know, 
for  application  purposes,  the  extent  to  which  aspect  variability  is  attribut¬ 
able  to  natural  causes  such  as  magnetic  field  variations. 

(S)  The  next  problem  is  the  determination  of  the  semi  thickness  A.  If  the 
height  h  ,  of  reflection  is  known,  the  easiest  way  to  determine  this  is  from 
plots  such  as  figures  9.5,  9.6,  9.7,  and  9.8.  However,  the  height  of  reflee 
tion  h  depends  on  the  angle'  of  incidence,  and  therefore  should  he  determined 
by  ray  tracing  for  each  individua1  profile  if  a  precise  h^  is  to  be  known. 

This  may  present  some  difficulty  as  it  is  not  clear  that  horizontal  stratifi¬ 
cation  can  be  assumed.  The  quantity  X  can  also  be  determined  in  the  course 
of  this  comparison. 

(S)  It  seems  that  plasma-line  observations,  when  available,  might  be  a  use 
ful  diagnostic  for  the  variation  of  with  latitude  and  also  for  determining 
the  north-south  offset  x^.  Some  attempt  has  been  made  to  use  the  published 
material  on  White  Sands  monostatic  observations  to  compare  with  the  predictions 
of  the  previous  section.  However,  there  is  a  dearth  of  accurately-calibrated 
A-scope  or  elevation  scans.  To  make  meaningful  comparisons,  these  scans 
should  be  accompanied  by  values  of  and  by  accurately  calibrated  ranges. 

Given  a  good  selection  of  such  scans  (including  times  when  only  a  very  small 
return  was  seen)  will  enable  the  parameters  of  the  model  to  be  evaluated  with 
confidence . 

9.5  Bistatic  Observations  (U) 

(S)  Attention  needs  to  be  given  to  the  question  of  evaluating  the  cross  sec¬ 
tion  for  bistatic  propagation  between  two  sites  at  different  azimuths  from  the 
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he;., ml  region.  The  straightforward  Bookf  r  formula  does  not  apply  under  these 
conditions.  It  is  presumably  a  simple  wavelength  equivalence  in  that  the 
spatial  wavelength  producing  bistatic  scatter  is  i/2  cos  (<f>/2j,  where  >1  is 
the  azimuthal  angle  of  the  scatter  around  the  magnetic-field  axis.  Therefore, 
the  cross  sections  appropriate  to  such  bistatic  scatter  would  be  appropriate 
from  monostatic  scatter  at  a  frequency  cos  0/2)  times  the  frequency  used. 

The  details  of  cnis  relationship  should  be  ar  integral  part  of  the  scattering 
model  . 
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10.  EVALUATION  OF  SMALL -STRUCTURE  MODEL  PARAMETERS  (IJ) 

(S)  As  an  example  of  the  tvpe  of  analysis  discussed  in  Section  9.4,  the  pre¬ 
dictions  of  the  simple  model  were  compared  with  data  provided  by  Dr.  P.  Fialer 
of  SRI.  Using  a  mean  latitude  of  34.4  deg  for  the  SRI  site,  and  assuming 
monostatic  operation,  the  curves  shown  in  Figures  10.1  through  10.4  were  gen¬ 
erated  showing  the  relative  cross  section  per  unit  area  for  heater  reflection 
heights  of  230,  23S,  240  ,  and  245  k®..  The  optimum  height  of  reflection  is 
about  225  km  for  a  magnetic  field  model  with  a  dip  of  68.2  deg  over  Platte- 
villc,  and  all  the  heights  calculated  are  above  that  altitude.  Therefore, 
the  observed  returns  are  towards  the  northward  side  of  the  heated  region,  re¬ 
sulting  in  a  slant  range  for  the  greatest  cross  section  that  increases  with 
altitude  of  heater  reflection.  Curves  are  shown  for  a  range  of  values  for 
the  thickness  2A  of  the  heated  region. 

(S)  Two  techniques  were  used  to  deduce  the  quantity  A.  In  the  first  of  these, 

the  range  extent  of  the  VHF  scattered  signal  for  a  given  height  of  reflection 

was  matched  against  the  appropriate  curve  from  Figures  10.1-10.4.  Lacking 

quantitative  data,  the  10-dB  points  on  the  cross  section  curves  were  assumed 

to  delimit  the  range  extent  of  the  observed  echo,  In  fact,  a  detailed  com- 

/ 

parison  with  the  A-scope  curve  would  give  much  more  meaningful  results. 

(S)  One  example  presented  by  Dr.  Fialer  is  as  follows: 

Time  (GMT)  Range  Spread  Relative  Range  Heater  Reflection 

of  Peak  Cross  Height 

Seccion 

0314  May  12  131  km  0  235  km 

0338  104  20  km  246 
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RELATIVE  CROSS  SECTION 


ELEVATION  (deg) 

17  16  15 


SLANT  RANGE  (km) 

Figure  10.1 

Relative  cross  section  vs  slant  range  or  elevation  showing 
the  effect  of  2A  for  hr  =  230  (U) 
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17  16  15 
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Figure  10.3 

Relative  cross  section  vs  slant  range  or  elevation  showing 
the  effect  of  2A  for  h  =  240  (U) 
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Figure  10.4 

Relative  cross  section  vs  slant  range  or  elevation  showing 

the  effect  of  2A  for  h  =  245  (U) 
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Using  Figures  10.2  and  10.4,  the  range  spreads  for  these  two  times  correspond 
to  values  for  2A  of  35  and  20  km,  respectively. 

lie  :  irameter  2A  was  then  redetermined  from  the  change  in  range  of  the 
maximum  cross  section.  Figure  10.5  shows  how  this  range  depends  on  heater 
reflection  height  h  for  =  0  (i.e.,  disturbed  region  centered  over  Platte- 
v i II  and  with  50  km  for  X  (the  e-folding  horizontal  radius  of  the  heated 
region).  Interpolation  from  this  curve  shows  that  a  20-km  increase  in  slant 
range  would  accompany  an  11-km  increase  in  heater  reflection  heights  for  a 
value  of  26  km  for  2A.  It  is  felt  that  this  result  is  in  reasonable  agreement 
with  that  deduced  by  the  first  method. 

(S)  A  note  of  caution  needs  to  be  sounded  Concerning  the  assumptions  made  in 
this  model.  The  rapidity  with  which  the  range  of  maximum  cross  section  varies 
with  heater  reflection  height  depends  on  the  slope  of  the  heater  antenna  pat¬ 
tern  (on  its  northern  edge,  in  this  case).  This  in  turn  depends  on  the  as¬ 
sumed  value  for  X.  As  an  example,  Figure  10.6  (for  X  =  100  km)  shows  sub¬ 
stantially  increased  slopes  from  the  corresponding  curve  of  Figure  10.5.  The 
assumption  of  a  1  ?ater  volume  centered  north  of  Platteville  (x^  =  30)  is  il¬ 
lustrated  in  Figure  10.7.  Its  principal  effect  is  to  displace  the  curves  of 
Figure  10.5  to  greater  ranges  and  greater  reflection  heights  without  aff  cting 
their  slopes . 

(S)  No  doubt,  a  detailed  comparison  of  these  curves  with  experimental  data 
will  permit  the  evaluation  of  all  three  of  the  parameters  A,  X,  x^.  However, 
present  indications  are  that  the  scattering  region  at  VHF  may  be  only  20  km 
thick,  implying  critic?!  geometry  for  any  bistatic  experiment . 
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SLANT  RANGE  (km) 


ELEVATION  (deg) 


SLANT  RANGE  (km) 


ELEVATION  (deg) 


Figure  10.7 

Slant  range  or  elevation  for  maximum  cross  section  vs_  h 
showing  the  effect  of  2A  for  x«  =  30  and  X  =  50  (U) 
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11.  LARGER-STRUCTURE  MODEL  (U) 


(S)  The  second  part  of  this  model  involves  the  interpretation  of  satellite 
scintillation  data  which  are  affected  by  artificial  spread  F  (ASF)  at  heights 
above  180  kin.  Because  these  irregularities  appear  to  extend  throughout  thr 

F  region,  the  height  p-.tent  of  the  irregularities  is  larger  than  the  horizontal 

» 

extent  of  the  disturbed  region.  Therefore,  the  disturbed  region  must  be  con¬ 
sidered  as  a  cylinder  tilted  at  an  angle  6  to  the  horizontal.  At  a  given  al¬ 
titude,  its  latitude  variability  is  as  shown  in  Equation  9.13;  however,  its 
height  variation  seems  most  rerdily  interpreted  as  being  due  to  a  constant 
percentage  fluctuation  in  electron  density  at  all  altitudes.  The  rms  phase 
fluctuation  y  is  given  in  Section  S  as 


5/2 


tH)' 


D  sec  i  cosec  /  (N/N  )^  dh 


(5.4) 


on  the  assumption  that  the  fractional  rms  electron-density  fluctuation  is  in¬ 
dependent  of  altitude.  The  rms  amplitude  fluctuation  e  is  given  by 


e  =  y//T  .  (5.5) 

The  rms  fractional  fluctuation  in  electron  density  was  found  in  Section  5  to 
be  about  6.9  x  10 

(S)  A  detailed  comparison  of  the  rms  amplitude  fluctuation  of  Figure  5.5 
through  the  satellite  pass  with  the  satellite  orbital  parameters  suggests  the 
form  shown  in  Figure  11.1  for  the  intersection  of  the  satellite  line  of  sight 
with  the  e-folding  contour  of  the  fractional  electron-density  fluctuation. 

This  corresponds  to  a  cylinder  of  radius  50  km  whose  axis  passes  through 
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the  300 -km  level  at  30  km  north  of  Platteville,  and  seems  a  useful  nominal 
mode  1  . 

(S)  Accordingly,  it  is  recommended  that  liquations  (5.4)  and  (5.5)  should 
be  used  to  determine  the  amplitude  fluctuation  expected,  using  the  following 
form  for  the  fractional  electron-density  fluctuation: 

(AN/N)2  =  (AN/N)2  exp  {[-(x-x0)2  -  y2 J/X2}  (11.1) 

where  x  and  y  are  distances  measured  in  the  northward  and  westward  direction, 
respectively,  from  the  magnetic  field  line  passing  above  Platteville  at  300-km 
altitude,  and  the  parameters  have  the  following  values: 

(AN/N)2  ^  6.9  x  10"6 
xQ  =  30  km 
X  =  50  km 

(S)  Fluctuations  in  angle  of  arrival  of  the  transmitted  wave  may  be  expected 
to  have  a  magnitude  of  about 


Ay/2/2  tt  D 


(11.2) 
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12.  GEOMETRY  AND  CROSS  SECTIONS  FOR  BISTATIC  SCATTER 


12.1  The  Bistatic-Geometry  Problem 

Consider  a  set  of  magnetic  coordinates  such  that  the  latitude  and  longi¬ 
tude  of  P  are  X^  and  0,  respectively.  Consider  propagation  from  a  point  1, 
on  the  earth's  surface,  to  a  point  0  at  height  hQ,  and  thence  to  a  point  2  on 
the  earth's  surface.  Let  6^  be  the  magnetic  dip  at  the  point  0.  Let  the  lati¬ 
tudes  and  longitudes  at  points  0,  1,  and  2  be  denoted  by  X  and  <p  with  the  ap¬ 
propriate  suffix. 

Consider  now  a  set  of  Cartesian  coordinates  with  the  z  axis  at  the  N 
magnetic  pole,  and  the  (x,  z)  plane  tne  meridian  through  P.  The  Cartesian 
coordinates  of  the  points  are: 


x^  =  R^  cos  Xj  cos  <J>^ 


y^  =  Rq  cos  X ^  sin  <J>j 


=  Rq  sin  X^ 


x 2  =  R0  cos  X2  C°5  *2 


y2  =  Rq  cos  \2  sin  4>2 


z2  -  Rq  sin  x2 


x0  ■  (Ro  *  ho)  cos  xo  cos  *0 


y0  ■  tRo  *  V  cos  xo  sln  ♦{ 


zo  ■  <Ro  *  ho>  sln  xo 
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where  R  is  the  radius  of  the  earth. 

The  magnetic-field  direction  at  the  point  0  is  defined  by  the  unit  vector 


B  =  [s in (6q  +  Aq)  ,  0,  -cos  (6Q  +  A(,)] 


The  vector  7^  joining  the  points  1  and  0  has  the  value 


10 


3  ix0  -  v  y0  -  yr  zo  -  zd 


Now  if  0^  is  the  angle  between  and  B  at  the  poin  0, 


Rj  cos  =  r1Q  •  B 


=  (xQ  -  Xj)  sin(6Q  +  A0)  -  (zQ  -  cos(6Q  +  AQ)) 


\  here  =  |'r^(J  |  =  slant  range  of  0  from  1. 

Similarly,  the  angle  0^  between  anc*  B  is 


cos  62  =  (xQ  -  x2  sin (6^  +  AQ)  -  (zQ  -  z^)  cos(6Q  +  AQ) 


For  specularity, 


0j  =  180°  -  02 


cos  0^  =  -cos  @2  • 


Therefore, 


x0  "  X1  Z0  Z1 

ir —  sln(6o  *  V  -  — ■ — cos  wo  *  V 


'1 


*1 
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x0  -  x2 


Z0  "  Z2 


+  R2  sin  (6o  +  A0)  -  R2  ct)S(6o  4  A0)  =  0 


x0  sin(S0  .  X0) 


1  1 

1  1 

R1  R2 

-  20  cos(50  *  V 

*r^ 

*  si"C«0  *  XQ) 


X1  X2 

_  .  _  |  +  cos(6()  +  y 


ZJ_  ZJ_ 

R1  +  RZ 


Substituting  for  xQ  and  zQ ,  the  first  line  becomes 


(R0  +  V 


J.  +  j_ 

R1  R2 

. 

cos  XQ  cos  (pQ  sin(60  +  ,.Q)  -  sin  cos(6Q  +  AQ) 


and  since  <  1  degree,  cos  <J>0  ~  1  and  this  becomes 


*  V 

'  \ 

1  1 

*  ,  , 

“1 

z 

2 

R1  *  R2 

.  4 

sin  S0  =  cos(60*  X0) 

*r 

■b) 

'  xi 
Ri 

-  sin(S0  *  V 

from  which  (R^  +  h^)  may  be  found.  The  difficulty  of  not  knowing  and  R„ 
disappears,  incidentally,  if  =  R^ .  Then  the  equation  becomes: 


2(R0  4  hQ)  sin  6Q  =  z2  +  z2)  cos(6Q  +  AQ) 


-  (Xj  +  x2)  sin  (6^  +  XQ) 
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The  suggested  procedure  is  as  follows: 

(1)  Identify  points  1  and  2. 

(2)  For  a  point  300  km  over  P,  find  and  R^. 

(3)  Solve  the  above  equations  for  h^  for  P  as  point  0.  If  h^  is  not 

about  300  km,  points  1  and  2  are  wrongly  located  for  successful 
bistatic  propagation.  Relocate  them. 

(4)  Find  h^  for  a  given  point  0;  recalculate  R^  and  R^  and  repeat. 

(5)  Choose  three  points  0  arranged  around  P,  and  fit  an  hp  surface  as 

outlined  In  the  following  section. 

12.2  Suggested  Procedux’e  for  Bistatic  Calculation  of  h^  (Orthogonality) 
Surface 


Refer  to  Figure  12.1.  Select  three  points,  A,  B,  and  C,  at  the  corners 
of  the  equilateral  triangle  centered  on  P’atteville  (P)  with  its  apex  to  the 


north,  and  with  sides  of  70  km. 
mula  (for  point  A): 


For  each  point, 


calculate  h 


Q 


using  the  for- 


hA„  R2  sin(X1  -  Aa.  -  6AJ  R:  sin(X2  -  A^  -  6AJ 

R0  R1  +  R2  sin  R1  +  R2  sin  V 


where  suffixes  1,  2,  and  A'  refer  to  the  transmitting  and  receiving  sites, 
and  to  the  point  on  the  h^  surface  above  A,  respectively;  R^  is  the  earth  radius 
R^  and  R0  are  the  ranges  from  1  and  2  to  k'\  6^  is  the  magnetic  dip  at  A'; 
and  X  is  latitude  in  magnetic  coordinates. 

Fit  a  plane  through  the  points  A',  B',  and  C "  and  take  this  as  defining 
the  hp  surface  in  the  region  of  interest. 

12.3  Bistatic-Cross-Section  Calculation 

Adopt  the  coordinate  system  shown  in  Figure  12.2. 
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Figure  12.1 


Orthogonality-surface  geometry 
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Then  the  direction  cosines  (1^,  np  ard  (1^,  r.^)  for  the  two  rays  1 

and  2  are : 

1^  =  cos  Tj  1 2  ~  cos  V 2  cos  6 

irij  =0  m2  =  "C0S  **2  S'*'n  ® 

n  =  sin  n2  =  sin 

Let  the  correlation  function  be  given  by 


2  2  2 
2T  2T  2L 

P(x,y , z)  =  e 


and  the  spatial  spectrum  by 


2  2 

T  1 


T2  2 
I  m 


.2  2 
L  n 


P (l»m,n,) 


(2tt)3/^T2L  e 


i 

: 


1 


Booker  gives  the  scattering  cross  section  as 


a 


v 


2  .  2 
7t  sm  x 


P{k(l2 


-  1^,  k(m2 


nij)#  k(n2  -  np} 


where  x  is  the  angle  between  the  electric  vector  of  the  incident  wave  and  the 
direction  of  the  scattered  wave  2.  Putting 


) 


the  cross  section  is: 
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,,  ,3/2  2  .  2 

(2 it)  it  sin  x 


^ir"T2L  e-(l/2)k2T2[(l2-l1)2  +  (m2-mi)2]  -  (l/2)k“L:(n2-n1) 


Let  s  be  the  distance  measured  along  the  magnetic  field  from  the  specular  point. 
Then 


2 

n^  =  sin  4*  +  s  cos  Y/R^ 
n2  =  sin  41  -  s  cos“4'/R2 

where  R  and  R2  are  the  slant  ranges  for  rays  1  and  2,  and  S'  is  the  specular 
angle  between  the  wave  normals  and  a  plane  normal  to  the  field.  Putting 

2  -1  -1 

n2  ~  nl  =  ”  C0S  ^^1  +  ^2  ^ 

in  the  cross  section  and  integrating  with  respect  to  s,  the  cross  section  per 
unit  area  normal  to  the  field  becomes: 


°A  = 


t  ,3/2  2  ,  2 

o  ds  =  i£2 - 

V  ,  4 


AN1 

2 

,  NJ 

2J2.  r  i  i  ^  +  fm  -m.)  ] 


T2L  e.(l/2)k2T  [(l^ip  +  (Vml 


2224 

-(l/2)kVs  cos44'(1  +  R2  ) 

e  ds 


a 


A 


an' 

2 

j  4  •  fc  rti  b 

4tt  sin  x  T 

l  nJ 

4  2-1  -1 

XN  k  cos  W(R1  +  R2  ) 


2  2 
-2kV 


cos2  V  sin2 (6/2) 
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Finally,  per  unit  area  of  horizontal  area  is  given  by: 


.  _  4TT4sir.2  x  c°s  <5 

AH  A„,4k  cos2'F(R1'1  +  R2' 


AN 

N 


,2„,2 


N 


2  t2£  - 2K  T  cos  V  sin  (6/2) 

/ 


For  a  medium  containing  more  than  one  scale  T,  the  terms  from 


AN 
.  N 


2 


onward  are  replaced  by: 


!— I2  V  2  e  -2k2'r12cos2'fsin2 (0/2)  an)  2  2  -2k 2T Aos^sin2  (6/2) 

l  N  J  T  1 1  e  1  T2  e  2 


where  6  =  magnetic  dip. 

12.4  Carrying  Qut  the  Integration 

Consider  a  coordinate  system  (x,y,z)  centered  at  Platteville.  A  suitable 
grid  of  prints  in  the  (x>y)  plane  is  selected,  and  the  geometrical  specular 
height  z  =  hq  is  calculated.  Then 


l~]  - 

[an 

l  N  J  ' 

{  N 

2  e_(x  ■  y2/x2  -  y2/*2  -C2  -  hQ)2M2 
0 


6  =  magnetic  dip  at  from  field  model,  and  Rj,R2»X*'fj6  are  calculated  at 

(x,y,hg).  Illumination  due  to  antenna  polar  diagrams  is  calculated.  The  pro¬ 
duct  is  integrated  over  the  (x,y)  plane. 

UNCLASSIFIED 

93 


SECRET 


13.  CONCLUSIONS  (U) 

13.1  Large-Structure  Model  (U) 

(S)  Observations  of  signals  from  geostationary  ami  orbiting  satellites  at  fre¬ 
quencies  from  30  to  400  MHz  have  established  the  general  behavior  of  the  large- 
scale  field-aligned  structure  (ASF)  that  produces  scintillations  of  those  signals. 
Studies  have  been  limited  to  full-power  ordinary-mode  heating,  with  the  heater 
being  directed  overhead,  at  frequencies  below  the  FT-layer  critical.  The  conclu¬ 
sions  are  as  follows: 

(S)  1.  The  field-aligned  structure  gives  strong  VHF  and  UHF  scintillations  up 

to  400  MHz. 

(S)  2.  The  scintillation  index  maximizes  when  the  direction  of  propagation 

is  close  to  the  downfield  direction.  Therefore,  for  example,  heating 
the  ionosphere  upfield  from  a  radar  will  produce  a  "blind  spot"  of 
about  10  deg  radius  centered  on  the  magnetic  zenith  at  150  MHz  frequency. 
(S)  3.  The  fading  depth  of  the  scintillation  varies  inversely  with  frequency, 

as  would  be  expected  from  the  integrated-phase-path  theory  of  satellite 
scintillations . 

(S)  4.  The  transverse  scale  of  the  large-scale  structure  has  a  half-correlation 

distance  of  about  100  m,  and  convects  in  an  east-west  direction  with 
the  velocity  of  fehe  ambient  plasma  (about  20-50  m  sec  1). 

(S)  5,  The  large-scale  structure  extends  from  200  to  4S0  km  altitude,  and  is 

usually  confined  within  a  cylinder  of  about  50  km  radius,  with  its 
axis  the  magnetic  field  line  passing  through  the  reflection  point  directly 
over  the  heater. 

(S)  6.  Occasionally,  unusual  shapes  for  the  irregular  region  are  seen,  probably 

arising  from  a  doughnut-shaped  disturbed  region  due  to  bum-through 
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>t  the  central  core  of  the  ionosphere  above  the  heater. 

(S)  Thc  f°llowing  additional  studies  of  large-scale  structure  are  important: 

(S)  1 .  Measurements  of  the  height  distribution  of  the  electron- density  fluctu¬ 

ation  associated  with  the  large-scale  irregularities.  Heights  inferred 
trom  orbital-satellite  passes  seem  to  point  to  field-aligned  irregular¬ 
ities  that  are  distributed  throughout  the  bottomside  and  topside,  with 
a  statistically  constant  fractional  fluctuation  of  electron  density. 
However,  this  has  not  been  experimentally  verified.  A  relatively  simple 
experiment  to  do  this  will  be  to  observe  the  rise  or  set  of  an  orbiting 
satellite  behind  the  disturbed  region,  as  observed  from  a  site  far  to 
the  south  (e.g.  Albuquerque).  Since  the  line  of  sight  to  the  satellite 
is  nearly  orthogonal  to  the  magnetic  field,  the  altitude  intercepted 
within  the  disturbed  cylinder  is  much  more  nearly  constant  than  when 

looking  at  an  orbital  satellite  from  a  location  to  the  north  of  the 
disturbed  region. 

(S)  2.  Detailed  studies  of  the  spatial  and  temporal  characteristics  of  the 

correlation  function  of  the  amplitude  for  an  orbiting  satellite.  Thus, 
it  will  be  possible  to  determine  whether  finer-scale  structure  is  pre¬ 
sent  on  the  scintillations.  Adequate  signal -processing  instrumentation 
is  available  to  make  these  studies  as  time  permits. 

(S)  3.  Studies  of  the  aspect  dependence  of  the  scintillating  signal  close  to 

the  down  field  direction.  Scattering  theory  predicts  the  way  the  ampli¬ 
tude  and  phase  fluctuation  should  vary  with  angle  from  the  magnetic 
zenith;  but  because  of  the  possible  significance  of  radar  applications, 
this  aspect  of  the  theory  needs  verification. 

(S)  4.  The  horizontal  morphology.  The  factors  controlling  the  shape  of  the 

disturbed  region  are  not  wholly  understood  at  this  time,  and  additional 
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morphological  studies,  perhaps  with  other  orbiting  satellites,  are 
appropriate. 

(S)  5.  Yield  studies.  Most  studies  up  to  now  have  concentrated  on  formulating 

a  scattering  model  appropriate  to  full  power  ordinary  polarization. 
Attention  needs  tc  be  given  to  the  effect  of  reducing  the  power,  opti¬ 
mizing  the  heater  frequency,  steering  the  heater  beam  and  trying  other 
forms  of  heater-signal  modulation. 

(S)  6.  Scintillation  experiments  at  ether  locations.  All  results  thus  far 

have  been  restricted  to  the  Platteville  heater  installation.  No  ASF 
has  been  demonstrated  at  other  locations.  Therefore,  it  seems  neces¬ 
sary  to  try  a  satellite  scintillation  experiment  at  Arecibo  (when  that 
heater  is  in  operation)  to  see  if  there  is  any  constant  difference  in 
the  large-scale  structure  between  Arecibo  and  Platteville.  In  addition, 
if  a  mobile  heater  facility  is  set  up  elsewhere  than  these  two  locations, 
satellite  scintillation  seems  to  be  the  appropriate  diagnostic  to  include. 


13.2  Smal  1- Structure  Model  (U) 

(S)  In  Sections  8,  9  and  10  of  this  report,  the  concept  has  been  developed  that 
the  small-scale  structure  responsible  for  VHF  and  UHF  backscatter  is  governed 
by  the  geometry  of  the  surface  of  orthogonality  (the  h^  surface)  appropriate  to 
the  location  of  the  radar.  The  intersection  of  this  (ideally  infinitely  thin) 
surface  with  the  disturbed  region  gives  the  location  of  the  backs cattering 
irregularities  for  that  particular  radar.  Using  a  simplified  magnetic  field 
model,  it  was  shown  that  the  h^  surface  slopes  down  toward  the  south,  and  that 
the  thickness  of  the  disturbed  region  can  therefore  be  deduced  from  the  range 
spread  of  a  backscattered  echo.  An  illustrative  example  was  given  to  show  that 
the  region  thickness  was  about  20  km.  A  methodology  has  been  indicated  for  the 
calculation  of  the  corresponding  geometry  and  cross  sections  for  bistatic  scatter. 
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(S)  The  following  important  tasks  remain  with  the  small-scale  scattering  model: 

(S)  1.  Disturbed  region  thickness  calculation.  The  same  type  of  calculation 

as  shown  in  Section  10  of  this  report  needs  to  be  carried  out  under  a 
wide  variety  of  conditions,  to  determine  the  limits  on  the  disturbed 
region  tnickness.  The  assumption  that  the  scattering  region  at  a  given 
range  is  infinitely  thin  needs  to  be  verified  by  an  interferometer  or 
other  similar  experiment. 

(S)  2.  Bistatic  geometry.  The  change  in  total  cross  section  with  scattering 

angle  indicated  in  Section  12  needs  to  be  verified  by  bistatic  tests. 

(S)  3.  'Held  studies.  Variations  of  the  small-scale  structure  cross  section 
and  morphology  with  frequency  and  reflection  height  and  heater  power 
need  to  be  evaluated. 
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